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Abstract 
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I. INTRODUCTION 



The contact interaction, or as sometimes called point or zero range interaction, plays 
an important role in nuclear physics. The short range nucleon-nucleon interaction is often 
simulated by a contact interaction in different models. An effective interaction of contact 
type (Skyrme-force) has been very frequently and successfully used in mean-field descriptions 
of heavy nuclei including the Hartree-Fock or Hartree-Fock-Bogoliubov theory (see e.g. Jl|). 
Recently contact interaction has been applied in three-body calculations of halo nuclei 
and in the study of the role of resonance states in pair correlations ||. Another example 
is the constituent-quark model for the application of the contact interaction. Several non- 
relativistic three-body calculations have calculated the spectra of the three quark system 
Despite the popularity of the contact force even its definition is problematic. The 
Hamiltonian of a two-body system interacting via a contact interaction is written as 



where /i is the reduced mass, A r is the three dimensional Laplacian, 5(r) is the Dirac 
delta-function and A is the "strength" of the contact interaction. The interpretation of the 
equation ([I]) is not trivial due to the singular nature of the Dirac delta-function. By using a 
complete set of states to diagonalize the above Hamiltonian one gets infinite binding energy 
for the system. To get a more meaningful result one usually truncates the configuration space 
@-|5||§ with some "reguralization" scheme. This truncation leads to a model dependence 
and results of different approaches can only be related and reproduced if one uses the same 
strategy to reguralize the interaction. Due to this approach it is widely believed that the 
contact interaction is meaningful only in a truncated space. 

In this paper we propose a different approach where the application of the contact in- 
teraction is based on mathematically rigorous theory |J. This proper contact interaction 
(PCI) is free from the above mentioned problems. The PCI gives finite energy and there is 
no dependence on basis truncations or extra parameters. This PCI has zero range and it 
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(i) 



acts non-trivially only in the I = partial wave. These two properties are what physicists 
expect from a contact interaction. 



The contact force was introduced in nuclear physics by Bethe and Peierls |L0j and Thomas 
I 1| in 1935. The mathematically correct treatment of the operator ([!]) dates back to the 



work of Berezin and Faddeev [12| in 1961. An English review of the early Russian works on 



two and three-body problems with PCI can be found in Ref. ||13|| . Mathematically rigorous 
treatment of the operator (|l|) is given by different approaches. One of them is the coupling 



constant re-normalization procedure of Berezin and Faddeev |12| . This technique was applied 
also in two dimensions ]14[]. Other approaches are the reguralization of the Green-function 
15,16], the dimensional reguralization |T7[ , or the theory of the Dirichlet forms We will 



use the self-adjoint extension theory of the contact force which is described in great detail 
in the monograph [[]]. Each of these approaches leads to the same spectrum of the two-body 
Hamiltonian with PCI. 

The mathematically correct formulation of the PCI in two-body problem is completely 
solved. It gives the two-body t-matrix in analytical form. This t-matrix can be used in 
subsequent three-body calculations based on the Faddeev equations. However, if we want to 
solve the three-body Schrodinger equation in coordinate space using some approximations 
then we need a local realization of the PCI. The aim of the paper is to give a few explicit 
local potential realizations of the PCI and apply it in three-body calculations. 

The plan of the paper is the following. In section II the two-body spectrum of the PCI 
is discussed. In part III the general theory of the local potential realization of the PCI is 
reviewed and a few explicit realizations are given. In sections II and III. A those theorems 
of the monograph || are reviewed which are relevant for our purpose. In section IV the 
numerical solution of the two-body problem with PCI are shown. Finally, at the section V 
we present three-body calculation with PCI for the halo nuclei 11 Li. 



II. CONTACT INTERACTION IN TWO-BODY PROBLEM 



Our review of the mathematically correct treatment of the contact interaction is based on 
the monograph of Albeverio et al. |§. The self-adjoint extension approach can be originated 
from the intuitive picture of the Dirac delta-function. The "function" 5(r) is zero if r is 
different from the zero vector 0. This means that the "operator" H$ is "identical" with the 
kinetic energy — |^A r if r 7^ 0. In mathematical terms this means that the operator — |^A r 
acts only on the functions which belong to the function space C^°(TZ 3 \ 0). This restricted 
Laplacian is not self-adjoint but its symmetric closure has many self-adjoint extensions. 
The self-adjoint extensions are characterized by a real number a (—00 < a < 00). This 
self-adjoint extension, H aj s, will give the proper representation of the contact interaction. 

The spectrum of the operator H a ^ is the following. The continuous spectrum is [0, 00) 
for each a. If a < then the I = partial wave component of H a ^ has one bound state 
with energy 

E a , s = - (Ana) 2 (2) 

The Hamiltonian H a j has no bound states for higher partial waves. If the parameter of the 
self-adjoint extension is non-negative then there is no bound state in any partial wave but 
there is one virtual state in the partial wave / = with energy given by Eq. (Q). 

The continuous spectrum may be characterized by the phase shift 8" ,5 (k), where k denotes 
the wave number k = (2fiE /Ti 2 ) 1 / 2 and E > 0. In the partial wave / = the phase shift 
reads 

kcot(6o' S (k)) = — , (3) 

<V<5 

where 

*•* = -iL- (4) 

The higher partial wave phase shifts are identically zero. The S-matrix in the partial wave 
/ = has the following form 



u v ' Ana -ik 

The full characterization of a Hamiltonian is given if we know its Green-function. The 
Green-function of PCI in coordinate space can be given in closed analytical form |J. 

The effective range expansion of an ordinary local spherical potential is of the form 
k 2l+l cot5i(k) = — 1/ai + r\k 2 : /2 + where a\ is the scattering length and r\ is the effective 
range. The comparison of Eq. (|3]) with this expansion shows that the operator H a g contains 
an interaction with zero range. The self-adjoint extension parameter not only determines 
the binding energy (see Eq. (fj) if a < 0)) but it is also connected to the scattering length 
by Eq. (|). 

At first it may seem strange that continuously many PCI exist. One can understand this 
using the following arguments (see e.g chapter 14.1 in Ref. [|L9[0. If a spherical potential 
vanishes for r > r$ then the solution of the radial Scrodinger equation is determined by 
assigning a given value to the logarithmic derivative of the wave function at r$. This bound- 
ary condition depends on the energy. If the range of the potential Tq tends to zero (contact 
force) then the boundary condition parameter becomes independent of the energy and in 
this way we can assign any value to it. Thus continuously many zero range interactions can 
be defined. This approach is an another way to define the PCI. 

To sum up, the two-body problem with Dirac delta-potential can be solved analytically 
in a mathematically rigorous way. Infinite binding energy due to the contact force does not 
occur. The mathematically correct treatment of the contact force always leads to a finite 
and well defined results. 



III. LOCAL POTENTIAL REALIZATION OF THE CONTACT FORCE 

The abstract PCI can be realized by an appropriate limit of local potentials. We have 
to construct this local representation in such a way that it corresponds to the self-adjoint 
extension theory. Naively taking a potential which resembles to a Dirac 5 function will not 
work. 



A. Realization of the contact interaction by local potential 



The recipe for local realization of the contact interaction is given in the monograph |9| 
(see theorem 1.2.5). Let us take any local potential V(r) which is in the Rollnik class PJT9|] 
and assume that / dr(l + r) |V(r)| is finite. We define a Hamiltonian family by 



H 5 (e) = ~A r + ^.V(r/e), (6) 



where e > and A(e) is an arbitrary analytic function of e such that A(0) = 1. Let us assume 
further that the potential V(r) has a non-degenerate zero energy resonance, i.e. 

h 2 

- — A r $ + V{r)<& = 0, (7) 

for some $. Because $ is a resonance state it is not square integrable function. The 
normalization is defined by the following condition 

drV{r) |$(r)| 2 = -1. (8) 

Under these circumstances the Green-operator of the Hamiltonian family (^) converges to 
the Green-operator of H a ,s in norm resolvent sense if e — * 0. The self-adjoint extension 
parameter is given by 

a =_A'(0)£-( f drV{v)®{v)\ . (9) 
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The theorem 1.2.5 of Ref. includes other cases. Here we have quoted only the part that 
we will use. 

B. Realization by Poschl- Teller potential 

We will realize the result of the previous section using the Poschl- Teller potential. It is 
known that in this case the radial Scrodinger equation 

" ^» " ^ S{S u!L <r) = Eu(r), (10) 
2/i 2/itr cosh (r/b) 2 



can be solved analytically P0 |. The energy eigenvalues are given by 

E n = -(s-l-2n) 2 ^, (11) 

where n = 0,2,4. . ., and s is an arbitrary complex number. If we choose the strength so 
that s = 2u + 1, where v — 0, 2, 4, . . ., then we will have solutions with zero energy for the 
state n = v. The corresponding wave functions are 

Mr) cc Cl!Xx (tanh(r/6)) , (12) 

1/2 

where C 2 l + \ is the Gegenbauer polynomial. It is easy to see that these solutions are not 
square integrable so we have potentials which have zero energy resonance states. For sim- 
plicity we will consider the case v = 0. The <3> of Eq. (0) with the normalization ([8]) now 
reads 

t(r) ^y /2 ^w>. (I3) 

h rv87r 



The self-adjoint extension parameter can be calculated using Eqs. @ and (0), 

Let us introduce the Hamiltonian family corresponding to Eq. (^), 

Hp T {e) = -— A r - 15 

2/i 6 2 e 2 2/icosh 2 (r/(6e)) 

We will call this operator the scaled Poschl- Teller realization of the PCI. According to the 
previous section the resolvent of ([15]) converges to the resolvent of H a ^ as e — > 0. We thus 
have a local realization of the PCI. 

In this example the convergence of the Hamiltonian family flT5] ) in the limit e — > can 
be easily shown using elementary calculations. Since the eigenvalue problem of the scaled 
Poschl- Teller operator ( |l~5f) can be written in the form (|10D the eigenvalues of ( |15|) read 

B „ (e , n) = - (Al_ _ (16) 



where 



s(e) 



-l + (l + 8A(e)) 



1/2' 



(17) 



and n — 0, 2, 4, Let us calculate Ept(€, n = 0) in the limit e — > 0. We get 



lim Epr(e, 0) = lim 



's(e)-l\ 2 ft 2 



V(0) 2 



(18) 



since in the limit e -> we have s(e) — > 1. Calculating the derivative of s(e) at e = we 
end up with 



i^M) = -1a<(o)^. 



(19) 



Thus we get the same energy as in the abstract theory of the PCI, i.e., if we substitute (|T 
into ©. 

The s-wave S-matrix of the scaled Poschl- Teller potential is of the following form 



2ikbe 



T(ikbe)T(\{2 + 8 (e) - iJfe&e))r(±(l - s(e) - ikbe)) 



(20) 



.27 r(-zA;foe)r(|(2 + s(e) + iAr6e))r(|(l - s(e) + zA;6e)) ' 
Using the identity T(l + z) = zY(z) of the Gamma function we can easily calculate the limit 



of the S-matrix 



i- nPTn \ i- l-s(e)+ikbe s(0)+ikb 

lim Sn \k,e) = lim — — — = — — . 

e-o v ; e-ol- s (e) - iA:6e -s'(0) - ifcfc 



(21) 



Taking into account Eq. ( 14] ) and s'(0) = 2A'(0)/3 we see that the limit agrees with the 
expression (Eh given by the general abstract theory of the PCI. 



C. Realization by exponential and square well potentials 

A local realization of the abstract result can be given not only with the Poschl- Teller po- 
tential but by any local potential satisfying the conditions of section III. A. Here we take first 
the exponential potential as an another example. The eigenenergies of the radial Scrodinger 
equation, 



ti 2 

— —u(r) — Vq exp(— r/b)u{r) = Eu(r), 



(22) 



can be obtained from the transcendental equation JT9 



where J v is the Bessel function. The condition of the zero energy resonance reads 



(23) 



* I 26^=0. 



(24) 



This equation defines the strength of the potential. The corresponding resonance wave 
function with the normalization given in section III. A is 



$(r) 



-J (6 exp(-r/26)), 



(25) 



where b = 2&(2 / u,V /^ 2 ) 1/2 . 

The scaled Hamiltonian family which converges to the PCI in the limit e — > is 

- — A r -y ^exp {-r be ). 

Note that the potential strength Vq is determined by Eq. fl24"|). We call the operator ([2"6]) 
the scaled exponential realization of the PCI. The self-adjoint extension parameter can be 
calculated using Eqs. (01) and fl2"5]) 



(26) 



a = -A'(0) 



,, n s ft 2 6Ji(6o 



drr exp(—r /b)J (b exp(— r/2b)) 



(27) 



2/x 4ttV 

The integral in (p7|) can be calculated numerically. 

The abstract PCI can also be realized by a properly chosen square well potential. This 
latter is defined in the usual way: V ®(r) = — Vq if r < R and V ®(r) = if r > R. The 
condition of the presence of a resonance with zero energy is that the strength of the square 
well is 

h 2 f(2A; + l)7r n 2 



Vn 



2fl 



R 



(28) 



where k = 0, ±1,±2, .... The properly normalized zero energy resonance wave function 
reads 



$(r) = i(2 7 rV i?)- 1 / 2 sin U^V r ) . (29) 



The self-adjoint extension parameter is given by 



a = - X '(0) 2 4^. (30) 
V J h 2 8tt K ' 



The Hamiltonian family which tends to the PCI is the following 



^-A r + ^Vf%), (31) 



where Vf Q (r) = —Vq if r < Re and V/ Q (r) = if r > Re. 

IV. NUMERICAL CALCULATION: TWO-BODY PROBLEM 

We have seen that the eigenvalue problem of the Poschl- Teller realization of the PCI can 
be solved exactly for any e. Since we will solve the three-body problem numerically, using the 



stochastic variational method p2| , it is essential to see how this numerical method performs 
for the solution of Eq. (|i~5f) or any other local realization of the abstract PCI. This numerical 
exercise is indeed a new test for the stochastic variational method since the potential of (|15D 
really resembles the Dirac delta-function. A few scaled Poschl- Teller potentials are shown 
in Fig. 1. for 6=1 fm. 

For the function A(e) we have chosen the following ansatz A(e) = 1 + Ce. With the 
parameter C we can change A'(0) which influences the value of the self-adjoint extension 
parameter. This latter uniquely determines the bound state energy and the scattering length 
of the PCI. 

For the two body problem we considered the bound neutron-proton system and the 
scattering states of the neutron-neutron system. For the deuteron case the exact eigenvalues 
and the numerically calculated ones using the scaled Poschl- Teller potential are compared in 
Table I for several values of the scaling parameter e. We have used h 2 / (2/x) = 41.47 MeV fm 2 



in all calculations. Table I demonstrates that the stochastic variational method is able to 
deal with such an extreme potential. Of course, numerically, we could not carry out the 
limit e — > 0. Using the reliable results of the stochastic variational method (e < 0.005) the 
numerically calculated values can be extrapolated to e = 0. From Table I we can conclude 
that the calculation with e = 0.005 is very satisfactory; the error is less than 10 eV. Having 
this experience, in three-body calculations we will accept the result of the calculation with 
scaling e = 0.005 (or the extrapolated value to e = 0) as the result of the PCI. 

The numerical results of the exponential and square well realization of the PCI are 
given in Table II. The parameters of the potentials were Vq =0.936831 MeV, b =8 fm and 
Vq =4.092925 MeV and R = 5 fm in the case of the exponential and square well potentials, 
respectively. The stochastic variational calculation performs very well for the exponential 
potential. It was not possible get better accuracy than four digits in the case of a square 
well potential. The convergence of the exponential and square well realization to the PCI is 
not as quick as in the case of the Poschl- Teller realization. 

The phase shift of the scaled Poschl- Teller potential 5(f T (/c, e) can be deduced from Eq. 
(pOl). The quantity k cot SQ T (k, e) is displayed in Fig. 2 as the function of the wave number 
for several values of the scaling parameter. We see that decreasing the scaling parameter 
the result converges to the phase shift of the PCI. The self-adjoint extension parameter of 
the PCI in this calculation is fixed to give the experimentally observed scattering length 
(—18.5 fm) of the dineutron. In order to get the effective range parameters of the scaled 
Poschl- Teller potential we make the expansion: 

k cot 5 PT (M = ik(S? T (k,e) + l)/(So PT (M - 1) « -l/a nn + r 2 k 2 /2. (32) 

At each scaling the scattering length and the effective range obtained in this way are shown 
in Tables I and III for the scaled Poschl- Teller potentials. 

To sum up, in numerical calculations we can only decrease the effective range to 0.01 fm 
but this is already enough to get the binding energy with an accuracy less than 10 keV. 



V. CONTACT INTERACTION IN THREE-BODY PROBLEM 



We have seen that the contact interaction is very tricky even for the two-body problem. 
Further complications can be expected when three particles interact. If each pair interaction 
is a contact force then the three-body problem can be solved exactly |]T2|,|T3J. The final result 
is, however, physically unsatisfactory due to the Thomas effect: there is an infinite number 
of bound states and the spectrum is not bounded from below |TTl , |2ll . We mention that in 
this three-body problem the self-adjoint extension type approach have to be used to have a 
meaningful three-body Hamiltonian. 

In the present paper we consider a three-body system where only one pair of the particles 
interacts through a contact force and the other two pair interactions are given by short range 
local potentials. The physical motivation of such a three-body system is to model nuclei 
with halos, like 6 He or n Li. The situation in this case is simpler. If we want to include the 
PCI in this type of three-body problem we can proceed in two different ways. 

For the solution of the three-body problem we may use the Faddeev theory. In this case 
one needs the Green-operator of each subsystem of the three particles. Since one pair of 
particles interact via PCI the Green-operator of this two body subsystem is known exactly 
and analytically. This property leads to a simplification of the Faddeev equations. It can be 
shown that this type of three-body problem is well defined and physically meaningful. The 
theoretical framework of this Faddeev type approach was developed in |23| . 



If we want to solve the three-body problem using the Scrodinger equation in coordinate 
space then we are ready to use the PCI since we have several local realizations of the PCI. 
Since we are interested in bound state properties of three-body systems we can use the 
Rayleigh-Ritz variational principle to calculate the binding energies. For the three body 
problem we have chosen a simple model of the light halo nucleus n Li which presently is in 
the focus of many theoretical and experimental investigations. 

Let us consider a simple model: an inert core and two neutrons. The three-body Hamil- 
tonian using Jacobi coordinates has the form 



h 2 h 2 

- ^A r - — A R + K,n(r) + V n>c {r, R). (33) 

Here r = ri — r2 is the relative coordinate of the two neutrons and R = |(ri + r 2 ) — r 3 is 
the vector which points from the third particle (core) to the center of mass of the dineutron. 
As in Refs. [|3],f|] we assume that the neutron-neutron interaction is a zero range contact 
force. However, we will not restrict the model space, as it was done in Refs. J3L|3J], in order 
to get finite energy but we will use the PCI. For the neutron-neutron interaction we will 
take the scaled Poschl- Teller and exponential realizations of the PCI. A shallow neutron 
core interaction is taken in order to avoid the problem of the Pauli projection. The explicit 
form of this potential is given in Ref. 0. 

In Table IV two calculations with the Poschl- Teller realization of the PCI are shown. In 
the first case the parameterization of the PCI corresponds to the experimental value (—18.5 
fm) of the scattering length of the free neutron-neutron system. The second calculation 
corresponds to the case when the parameter of the PCI is adjusted until the three-body 
energy reproduces the experimental value (S2 n — —0.295 ± 0.035 MeV) of the two-neutron 
separation energy of n Li p3j. This parameter of the PCI leads to a —7.32 fm scattering 
length for the dineutron. The mean square radius of the neutron-neutron and the core- 
dineutron system is also calculated in this case. 

From Table IV we can conclude that the effective neutron-neutron interaction is not 
consistent with the free neutron-neutron scattering. We had to reduce the scattering length 
of the dineutron system in order to reproduce the experimental value of the two neutron 
separation energy. Similar phenomena were observed in the calculations of Refs. where 
the traditional notion of the contact force was applied. The results of Table IV are very 
similar to the results of Ref. [|J where the standard method (restricted model space) was 
used to handle the contact interaction. Although the results are qualitatively similar there 
is a deep conceptual difference between the two calculations. In our calculations with PCI 
we indeed used a zero range force whereas in the works [0,|J, due to the restriction of the 
model space, the effective range is much larger than ours. In principle the PCI has zero 



range, in practice we can reach 0.01 fm. In the calculations of Ref. [[],|| the effective range 
is roughly 2 fm and this explain the quantitative difference of the two calculations. 

We have repeated the three-body calculation, replacing the Poschl- Teller realization by 
the exponential realization of the PCI. The results are shown in Table V. Of course, due to 
the different realization of the PCI, the energy corresponding to a given e are different in 
Table V and in the second column of Table IV. What is important is that the extrapolated 
values, which correspond to the PCI, are very close to each other in both realizations. 

Because we use the local realization of the PCI we have to carry out the limit e — > 
numerically in the three-body calculation and this requires great care. The stochastic 



variational method ||22|| , which proved to be very accurate in few body problems, is used to 
solve the bound state three-body problem. It is clear from Tables IV and V that convergence 
is reached in the order of a few tenth keV. This means that our proposed method for the 
correct handling of the contact interaction works for three-body problems. 



VI. SUMMARY 

We have reviewed the mathematically, rigorously defined proper contact interaction and 
shown that it has the properties of the standard notion of the contact interaction. This 
proper contact interaction is free from the problem (infinite binding energy) of the standard 
treatment of the contact interaction. Local potential realizations of the abstract proper con- 
tact interaction were given in terms of Poschl- Teller, exponential and square-well potentials. 
The calculation for 11 Li revealed that the proper contact interaction with local realization 
can be incorporated into three-body calculations based on the coordinate space Schrodinger 
equation provided this latter is solved accurately enough. 
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TABLES 

TABLE I. The binding energy in MeV of the deuteron in scaled Pdschl- Teller potential as 
the function of the scaling parameter. The exact result and the numerical one calculated by the 
stochastic variational method are shown. The binding energy of the proper contact interaction is 
2.222496 MeV this corresponds to the limit e — > 0. The effective range ro of the scaled Poschl- Teller 
potentials are also given. 



e 


r (fm) 


—Ept numerical 


—Ept exact 


0.500 


0.86 


2.065933 


2.065933 


0.100 


0.19 


2.188844 


2.188843 


0.050 


0.10 


2.205511 


2.205510 


0.010 


0.02 


2.219079 


2.219073 


0.005 


0.01 


2.220789 


2.220783 


extrapolated to 
e = 


0.00 


2.222502 


2.222496 



TABLE II. The binding energy in MeV of the deuteron in scaled exponential and square-well 
potentials as the function of the scaling parameter. The exact result and the numerical one calcu- 
lated by the stochastic variational method are shown. The binding energy of the proper contact 
interaction is 2.222496 MeV this corresponds to the limit e — > 0. 





e 


xponential 




square- well 




e 


numerical 




exact 


numerical 


exact 


0.500 


1.34894 




1.34894 


1.662 


1.663 


0.100 


1.93918 




1.93919 


2.077 


2.080 


0.050 


2.06796 




2.06797 


2.135 


2.149 


0.010 


2.18910 




2.18909 


2.196 


2.207 


0.005 


2.20562 




2.20562 


2.201 


2.215 


extrapolated to 












e = 


2.22250 




2.22250 


2.205 


2.222 


TABLE III. 


The effective range parameters of the dineutron system in scaled Poschl- Teller 


potential as the function of the scaling parameter. The effective range parameters of the proper 


contact interaction (scattering length 


arm = —18.5 fm and effective 


range ro = fm) corresponds 


to the limit e — > 


0. 










e 


0.500 


0.1 


0.05 


0.01 


0.005 


-a nn (fm) 


19.85 


18.38 


18.43 


18.49 


18.49 


r {fm) 


1.04 


0.20 


0.10 


0.02 


0.01 



TABLE IV. Two neutron separation energy of n Li in the three-body model as the function 
of the scaling parameter of the Poschl- Teller potential. The binding energy of the proper contact 
interaction can be obtained extrapolating to e = 0. The self adjoint extension parameter of the 
proper contact interaction corresponds to the given neutron-neutron scattering length a nn . The 
mean square neutron-neutron and core-dineutron radius are also given. 





a nT1 = —18.5 fm 


a nn = -7.32 fm 


€ 


S 2n (MeV) 


S 2n (MeV) 


rL ( fm2 ) 


r 2 cn (fm 2 ) 


0.500 


-0.5495 


-0.1590 


86.3 


113.0 


0.100 


-0.6949 


-0.2812 


58.0 


79.6 


0.050 


-0.7135 


-0.2975 


55.1 


76.0 


0.010 


-0.7268 


-0.3095 


54.3 


75.6 


0.005 


-0.7285 


-0.3112 


54.1 


75.3 


extrapolated to 










e = 


-0.7302 


-0.3130 


53.9 


74.9 



TABLE V. Two neutron separation energy of 11 Li in the three-body model as the function 
of the scaling parameter of the exponential potential. The binding energy of the proper contact 
interaction can be obtained extrapolating to e = 0. The self adjoint extension parameter of the 



proper contact interaction corresponds t( 


) the fr 


ee neutron-neutron scattering length —18.5 fm. 


e 




S 2n (MeV) 


0.100 




-0.2786 


0.050 




-0.4797 


0.010 




-0.6795 


0.005 




-0.7058 


0.0025 




-0.7188 


extrapolated to 






e = 




-0.7316 
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FIG. 1. Local potential realization of the proper contact interaction in the case of deuteron. 
Scaled Poscl-Teller potentials with different scaling parameters. Note the logarithmic scale on each 
axis. 
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FIG. 2. The phase shift of the scaled Poscl-Teller potential with different scaling parameters 
for the dineutron system. The parameter of the PCI corresponds to the experimental scattering 
length a nn = —18.5 fm. The quantity k cot(5o(k)) is displayed. For the proper contact interaction 
this is independent of the wave number k ( k cot(5o(A;)) = ffn -1 ) 



